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ABSTRACT: The swine leukocyte antigen (SLA)-DRA locus is noteworthy among other SLA class II loci for its limited variation and 
has not been investigated in depth. This study was investigated to detect polymorphisms of four exons of SLA-DRA gene and its 
association with piglet diarrhea in Landrace, Large White and Duroc pigs. No polymorphisms were detected in exon 3, while 2 SNPs 
(c.l78G>Aand c.211T>C), 2 SNPs (c.3093A>C and C.3I04OT) and 5 SNPs (c.4I67A>G, c.4I84A>G, c.4I94A>G, c.4246A>G and 
c.4293G>A) were detected in exon 1, exon 2 and exon 4 respectively, and 1 SNP (c.4081T>C) in intron 3. Statistical results showed that 
genotype had significant effect on piglet diarrhea, individuals with genotype BC had a higher diarrhea score when compared with the 
genotypes AA, AB, AC and CC. Futhermore, genotype AC had a higher diarrhea score than the genotype CC in exon 1 (p<0.05); 
diarrhea scores of genotype AA and BB were higher than those of genotypes AC and CC in exon 2 (p<0.05); individuals with genotype 
AA had a higher diarrhea score than individuals with genotype AB and BB in exon 4 (p<0.05). Fourteen common haplotypes were 
founded by haplotype constructing of all SNPs in the three exons, its association with piglet diarrhea appeared that Hap2, 5, 8, 10, and 
14 may be the susceptible haplotypes and Hap9 may be the resistant haplotype to piglet diarrhea. The genetic variations identified of the 
SLA-DRA gene may potentially be functional mutations related to piglet diarrhea. (Key Words: Swine Leukocyte Antigen, DRA Gene, 
Association, Genetic Polymorphism, Haplotype, Piglet Diarrhea) 



INTRODUCTION 

The major histocompatibility complex (MHC) class II 
molecules are tightly linked and highly polymorphic cell 
surface glycoproteins, which preferentially recognize 
foreign antigen to T cells and subsequently, initiate specific 
immune responses against infectious disease (Kedrick and 
Kim, 2000; Patch et al., 2011). In swine, swine leukocyte 
antigen (SLA) class II antigens are heterodimeric proteins 
that are comprised of non-covalently linked a (34 kDa) and 
P (29 kDa) chains, with the a chain coded by the SLA-DRA 
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and SLA-DQA genes and P chain coded by the SLA-DRBl 
and SLA-DQBl genes (Lunney et al., 2009). 

It is evidenced that some SLA alleles and haplotypes are 
associated with some immune-related diseases such as 
antibody responses to sheep erythrocytes and cell-mediated 
immune responses to tuberculin-purified protein (Mallard et 
al., 1989; Shinkai et al., 2012), PRRS (Lewis et al., 2007; 
Molina et al., 2008), and other swine immune response to 
various infectious agents (Renard and Vaiman, 1989; Chen 
et al., 2005; Ho et al., 2010a). To investigate the role of 
SLA genes in immune responsiveness, it is essential to have 
a better understanding of SLA polymorphism. 

The SLA-DRA gene, which encodes the a chain of DR 
antigen is closely linked to DRB. The full-length genome 
sequence of 4.6 kb, encodes for a polypeptide of 252 amino 
acid residues, which includes four exons and is considered 
to have limited allelic variation (Ho et al., 2009; Lunney et 
al., 2009). Therefore previous studies on genetic 
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polymorphism of SLA class 11 genes have primarily focused 
on polymorphism in DQB, DRB, and DQA genes 
(Luetkemeier et al., 2009; Ho et al., 2010a; Essler et al., 
2013), but rarely targeted the swine DRA gene for 
comprehensive analysis. Consequently, the extent of 
variation at the locus is not yet known. 

Due to large scale intensive production systems, pigs 
are exposed to various pathogens that cause infectious 
diseases and pathological symptoms (Lunney et al., 2009). 
Diarrhea in newborns and pre-weaning piglets, which 
causes increased mortaUty rates, poor growth and additional 
medical costs, is one of the most serious threats in the swine 
industry. It is estimated that 50% of piglet mortality is due 
to diarrhea deaths (Morris et al., 2002). Our preUminary 
studies indicated that piglet diarrhea during pre-weaning 
was linked to the exon 2 of SLA-DRA and SLA-DQA gene 
polymorphism (Yang et al., 2012; Yang et al., 2013). 

The objective of this study was to detect the genetic 
variation of the four exon regions of the SLA-DRA gene, 
and investigate the association with piglet diarrhea from 
Landrace, Large White and Duroc pigs using polymerase 
chain reaction- single strand conformation polymorphism 
(PCR-SSCP) and sequencing, to provide a theoretical 
foundation for further pig anti-disease breeding. 

MATERIALS AND METHODS 

Animals and genomic DNA extraction 

A total of 216 pure strain piglets from 24 sires were 
selected from Linze Xinghua Pig Breeding Farm (Zhangye, 
Gansu, China), including Landrace (n = 72), Large White (n 
= 72) and Duroc (n = 72). 

Experimental piglets were born in the same time period 
(within 20 days) and the experiment conducted under forage 
and feeding management conditions and disease control in 
accordance with the Guide for Chinese Feeding Standard of 
Swine and approved by the institute of Ministry of 
Agriculture Feed Industry Center, China Agricultural 
University (Beijing, China). All piglets were monitored two 
times a day during the entire suckling period (from birth to 
weaning, 0 to 28 days) and assigned a daily score based on 



visual analysis of symptoms traits as follows: 0: normal, 
soUd feces; 1: slight diarrhea, soft and loose feces; 2: 
moderate diarrhea, semi-liquid feces; 3: severe diarrhea, 
liquid and unformed feces (Kelly et al., 1990). 

Genomic DNA was isolated from swine ear tissues 
samples by standard phenol-chloroform extraction meth 
odology (Sambrook and Russell, 2001) and quantified 
by using a Nanodrop ND-1000 spectrophotometer (Ther 
mo Fisher Scientific, Wilmington, DE, USA). 

Polymerase chain reaction amplification and 
polymerase chain reaction-single strand conformation 

polymorphism analyses 

Four primer pairs (Table 1), were designed to amplify 
exons 1, 2, 3 and 4 of SLA-DRA gene according to the 
published SLA-DRA gene sequence in GenBank (accession 
No. AY303990). PCR analyses were performed in 25 ^L 
reactions containing of 100 ng genomic DNA, 2.5 ^L of 10 
xPCR buffer (containing Mg^"", 15 mM), 2.5 mM dNTP 
mixture, 10 pmol of each primer, and 5 U of Taq 
polymerase (TaKaRa, Dalian, China). The thermal profile 
consisted of 2 min at 94°C, followed by 35 cycles of 20 s at 
94°C, 30 s at annealing temperatures (Table 1) and 30 s at 
72°C, with a final extension of 10 min at 72°C. 

For SSCP analysis, 2 |iL PCR products were mixed 
with 7 nL of the denaturing solution (98% formamide, 
10 mM ethylene diamine tetraacetic acid (EDTA), 0.0 
25% bromophenol blue, and 0.025% xylene-cyanole), af 
ter denaturation at 98°C for 10 min, samples were coo 
led on ice immediately and loaded on polyacrylamide 
gel (170x170x1 mm) in IxTBE under optimal conditio 
ns. Gels were silver- stained according to the method of 
Byun et al. (2009). 

Cloning and sequencing 

PCR products representative of different SSCP patterns 
were purified using DNA Fragment Quick Recover Kit 
(TaKaRa, Dalian, China). The pmfied PCR products were 
Ugated into the pMD19-T vectors (TaKaRa, Dalian, China) 
according to the manufacturer's recommendations. The 
mixture was incubated at 4°C overnight, and transformed 



Table 1. Primer sequences, annealing temperature (Ta) and size of amplification products 


Locus 


Primer sequence (5'— >3') 


Ta (°C) 


Position 


Size (bp) 


Exon 1 


F: CTTTGCTTGTATTGC 
R: ACCTAACTACCCCTC 


56.8 


79-264 


186 


Exon 2 


F: CAGAGAATCACGTGATCAT 
R: ACAGGTACCATTGGTGTT 


55.0 


288-3137 


257 


Exon 3 


F: TGCTAAACAGGGAAGGCT 
R: ACAAAGGAGACTGAGGGATG 


56.8 


3531-3882 


352 


Exon 4 


F: TCCCGTAATACATCGTTC 

R: rrcc i'i icci roGC'ic Al' 


55.6 


4018-4374 


357 



Primer positions refer to SLA-DRA gene sequence in GenBank (Accession No.AY303990), exon 1, 2, 3, and 4 are located at 151-226, 2884-3129, 3573- 
3854, 4155-4309, respectively. 
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into Escherichia coli DH5a strains. Recombinant clones 
were detected by blue/white screening, 10 to 15 insert- 
positive clones for each transformation were incubated 
overnight at 37°C, in a vibrating rotary incubator (200 rpm). 
Verification of targeted inserts were screened using the 
PCR-SSCP approach as described by Zhou et al. (2008), 
and only those clones for which the SSCP patterns matched 
those of the corresponding genomic DNA were selected and 
subsequently sequenced by Sangon Biotech (Shanghai, 
China) Co., Ltd. 

Statistical analysis 

Nucleotide and amino acid sequences were analyzed 
using MEGA5 (Tamura et al., 2011) and new allele 
sequences were submitted to GenBank database 
(http://www.ncbi.nlm.nih.gov/nuccore/); gene heterozygosity 
{He) was calculated using POPGENE software (version 
3.2); polymorphism information content {PIC) was 
calculated by Botstein's methods (1980). Haplotype 
construction and linkage disequilibrium were analyzed by 
PHASE 2.0 (Stephens and Donnelly, 2003) and DNAsp 5.0 
(Hudson, 2001) software respectively. 

The effects of genotype on piglet diarrhea were 
analyzed using the general linear model procedure of SAS 
software (SAS, 1999) according to the following statistical 
model: 

^ ijklm 

^+Bi+SJ+Gk+eijkl 

Where Yyi^i is the observation value of piglet diarrhea 
score, |i is the overall population means, Bi is the fixed 
effect of breed, Sj is fixed effect of gender, Gk is the effect 
of genotype, and e is the random residual effect. 
Preliminary analysis also included the fixed interaction 
effects of breed, sex and genotype, however, these 
interaction effects were subsequently removed because they 
did not have a significant effect. 



RESULTS 

Genetic polymorphism of the SLA-DRA gene 

There was no polymorphism in exon 3 of the DRA locus. 
Six genotypes were detected in exon 1 and exon 2 (Figure 
lA and Figure IB), respectively. They were defined as AA, 

AB, BB, AC, CC, and BC. For exon 1, genotypes AA, BC, 

AC, and CC were detected in all three pig populations, 
genotype BB was only detected in Duroc. Three alleles 
were described as: A (A-T), B (G-C), and C (G-T) in this 
locus (Table 2). Allele C was the most common in Large 
White with a frequency of 0.55, whereas allele A was 
predominate in Landrace and Duroc with frequencies of 
0.5 1 and 0.44, respectively. 

For exon 2, genotypes AA, AB, BB, and BC were 
detected in all three pig breeds, but no AC or CC genotype 
were observed in the Duroc. Three alleles were described 
as: A (A-C), B (C-C), and C (C-T) at this locus (Table 3). 
Allele B was predominant in Large White and Landrace 
with frequencies of 0.51 and 0.61, respectively. The 
frequency of allele A (0.49) was identical with allele B in 
Duroc. 

Seven different genotypes (AA, AB, BB, BC, BC, BD, 
CE, and CF) were identified in exon 4 of the DRA locus 
(Figure IC). Seven genotypes were evident in the Large 
White and Landrace, but only genotypes AA, AB, and CE 
were detected in Duroc. Seven alleles were described as: A 
(T-A-A-A-G-G), B (T-G-A-A-G-G), C (C-A-A-A-G-G), D 
(T-A-A-A-A-G), E (T-A-A-G-G-A), and F (T-G-G-A-G-G) 
(Table 4). A was the most common allele in the three breeds 
with frequencies of 0.52, 0.54, and 0.90, respectively. 

Results of genetic indexes {He and PIC) indicated that 
all PIC values ranged from 0.4 to 0.59 (Tables 2, 3, and 4) 
with the exception of low polymorphism in Duroc in exon 4 
site {PIC = 0.19) (Table 4) indicating all of these loci 
exhibited abundant genetic diversity within the different pig 
strains according to the protocol of Vaiman et al. (1994). 




c 




BC AB AA BB CC AC AB CC BB BC AC BC AA CE BB AA AB CF BC BD 

Figure 1. PCR-SSCP analysis of PCR amplified products for exon 1, exon 2, and exon 4 of the SLA-DRA gene, respectively. Exon 3 with 
no polymorphism detected in the amplification region is not listed in this figure. (A) PCR-SSCP band patterns for exon 1 on 12% 
acrylamideibisacrylamide gel (39:1) after 10 min (250 V) and then 20 h (200 V) of the run at 4°C, (B) PCR-SSCP band patterns for exon 
2 on 10% acrylamideibisacrylamide gel (29:1) after 30 min (250 V) and then 20 h (200 V) of the run at 4°C, (C) PCR-SSCP band 
patterns exon 4 on 10% acrylamide:bisacrylamide gel (29:1) after 30 min (250 V) and then 20 h (200 V) of the run at 12°C. PCR-SSCP, 
polymerase chain reaction-single strand conformation polymorphism. 
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Table 2. Genotype and allele frequency, He and PIC of the SLA-DRA exon 1 in the three pig strains 







GenBank (No.AY303990) 5'^3' 




Frequency 




Exonl 




178 


211 


Large White 


Landrace 


Duroc 






^A 


fflCA 


72 


72 


72 




A A 


A 


T 


n 17 n ?■> 


n 1 ^ (()\ 

v.ij yy ) 


n 1 c n Q'l 




AR 


A/G 


T/C 


o.ui v^i } 


0 flfi C4~> 

U.UU ) 


\J.DJ } 




RR 










fl 1 fl C71 




BC 


G 


T/C 


0.04 (3) 


0.04 (3) 


0.11 (8) 




AC 


A/G 


T 


0.50 (36) 


0.71 (51) 


0.17 (12) 




CC 


G 


T 


0.28 (20) 


0.07 (5) 


0.10(7) 


Allele 


A 


A 


T 


0.42 


0.51 


0.44 




B 


G 


C 


0.03 


0.05 


0.33 




C 


G 


T 


0.55 


0.44 


0.24 


He 








0.51 


0.54 


0.67 


PIC 








0.40 


0.44 


0.57 



He, observed heterozygosity; PIC, polymorphism information content; SLA, swine leukocyte antigen. The following marker means are same. 
Numbers in brackets are population size tested that belong to the respective genotypes. Bases in the boxes are mutation bases in amplified regions. 



The PIC and He for Duroc pigs were lower than Large 
White and Landrace in exon 4 locus, which is not 
unexpected. However, it may be explained by the 
differences in genetic backgrounds as well as selection 
pressures between Duroc and Large White, Landrace. 

Nucleotide and amino acid variation 

The nucleotide characterizations were investigated by 
cloning and sequencing. According to references sequence 
AY303990 of SLA-DRA gene, two SNPs (c.l78G>A and 
c.211T>C) and two SNPs (c.3093A>C and c.3104C>T) 
were detected in exon 1 and exon 2, respectively (Tables 2 
and 3). Additionally in exon 4 locus, one SNP (c.4081T>C) 
in intron 3 and five SNPs (c.4167A>G, c.4184A>G, 
c.4194A>G, c.4246A>G, and c.4293G>A) in exon 4 were 
detected (Table 4). When compared with the existing SLA- 
DRA alleles in the Immuno Polymorphism Database (IPD)- 



MHC SLA sequence database (http://www.ebi.ac.uk/ipd/ 
mhc/sla/), mutations c.l78G>A, c.211T>C, c.4081T>C 
c.4184A>G and c.4293G>A were found for the first time in 
swine, which formed four new allele sequences and have 
been submitted to the GenBank database 
(http://www.ncbi.nlm.nih.gov/nuccore/) and assigned the 
accession number JX135565, JX135566, JX135567, and 
JX135568. 

Effects of genotypes and haplotypes on piglet diarrhea 

Results shows differences at SLA-DRA exon 1 and exon 
2 loci for genotype and piglet diarrhea (p<0.05) (Table 5). 
The least square means along with standard errors for piglet 
diarrhea score among different genotypes are given in Table 
6. 

In exon 1, the individuals with genotype BC showed a 
higher diarrhea score (1.77+0.28) than those with genotypes 



Table 3. Genotype and allele frequency, He and PIC of the SLA-DRA exon 2 in the three pig strains 







GenBank (No.AY303990) 5'^3' 




Frequency 




Exon2 




3093 


3104 


Large White 


Landrace 


Duroc 






gTG 


cg|c| 


72 


72 


72 


Genot5?pe 


AA 


A 


c 


0.13 (9) 


0.08 (6) 


0.21 (15) 




AB 


A/C 


c 


0.33 (24) 


0.33 (24) 


0.57 (41) 




BB 


C 


c 


0.26 (19) 


0.38 (27) 


0.18 (13) 




BC 


c 


C/T 


0.17 (12) 


0.14(10) 


0.04 (3) 




AC 


A/C 


C/T 


0.07 (5) 


0.06 (4) 


0(0) 




CC 


c 


T 


0.04 (3) 


0.01 (1) 


0(0) 


Allele 


A 


A 


C 


0.33 


0.28 


0.49 




B 


C 


C 


0.51 


0.61 


0.49 




C 


c 


T 


0.16 


0.11 


0.02 


He 








0.61 


0.54 


0.52 


PIC 








0.53 


0.47 


0.40 



He, observed heterozygosity; PIC, polymorphism information content; SLA, swine leukocyte antigen. 
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Table 4. Genotype and allele frequency, He and PIC of the SLA-DRA exon 4 in the three pig strains 









GenBank (No. AY303990) 5'- 


-^3' 






Frequency 




Exon4 




4081 


4167 


4184 


4194 


4246 


4293 


Large White 


Landrace 


Duroc 






T 


CgA 




a|a|c 


gg|a| 




72 


72 


72 




A A 


T 
1 


A 


A 


A 


u 


c 




0 4^ m 'I 


U.OJ VOJ^ 




AR 


T 
1 


A/G 


A 


A 


V T 


a 

u 




u.z..^ V^o_/ 


0 07 ^7^ 




RR 


T 
1 


G 


A 


A 




c 


u.u / ^ 


0 06 (A\ 


u (,u^ 




Rr 


TIC 


A/G 


A 


A 






n 1 n (i\ 


0 Oft rzl'i 


u tu; 




RFi 

lily 


T 
1 


A/G 


A 


A 


A/n 




n ni ni 

o.ui v^i } 


0 01 

u.ui v^i } 


0 CO'l 






TIC 


A 


A 


A/G 


\j 


n/A 


O.UO V_U/ 


u. lo V^''/ 






PF 
^1 


TIC 


A/G 


A/G 


A 


V. J 




0 04 (V\ 


0 04 C31 


0 ('O'l 


Allele 


A 


T 


A 


A 


A 


G 


G 


0.52 


0.54 


0.90 




B 


T 


G 


A 


A 


G 


G 


0.30 


0.20 


0.05 




C 


C 


A 


A 


A 


G 


G 


0.11 


0.14 


0.03 




D 


T 


A 


A 


A 


A 


G 


0.01 


0.01 


0 




E 


T 


A 


A 


G 


G 


A 


0.04 


0.09 


0.03 




F 


T 


G 


G 


A 


G 


G 


0.02 


0.02 


0 


He 
















0.63 


0.64 


0.19 


PIC 
















0.57 


0.59 


0.19 



He, observed heterozygosity; PIC, polymorphism information content; SLA, swine leukocyte antigen. 
SNP c.4081T>C was found in intron 3. 



Table 5. Effects of different factors on piglet diarrhea score 



Factors 




Exonl 






Exon2 






Exon4 




Breed 


Sex 


Genotype 


Breed 


Sex 


Genotype 


Breed 


Sex 


Genot5T3e 


F value 


0.65 


0.01 


3.21 


0.34 


0.02 


2.12 


0.65 


0.03 


1.58 


P 


0.52 


0.93 


0.008* 


0.71 


0.90 


0.04* 


0.52 


0.86 


0.16 



* Significant at 0.05 level. 



AA (0.69+0.18), AB (0.85+0.21), AC (0.97+0.11) and CC 
(0.52+0.19), and genotype AC showed a higher diarrhea 
score (0.97+0.11) than the genotype CC (0.52+0.19) 
(p<0.05). In exon 2, diarrhea score of genotypes AA 
(1.12+0.19) and BB (1.07+0.14) were significantly higher 
than those of AC (0.22+0.35) and CC (-0.40+0.53) 
(p<0.05). In exon 4, although the observed significance 
level for the F test was 0.1552, the genotypic linear 
contrasts suggested that diarrhea score was higher in the AA 
genotype (1.06+0.10) than AB (0.59+0.16) and BB 
(0.29+0.35) genotypes (p<0.05). The animals with genotype 



BD were excluded when association analysis was carried 
out because the sample size was less than three. 

Haplotype construction of all SNPs in the three exons 
and its association with piglet diarrhea are analyzed in Table 
7. A total of 54 haplotypes were found, 14 of which had 
frequencies greater than 2% and were defined as common 
haplotypes while the others were combined in this study 
(Breast Cancer Association Consortium, 2006). Diarrhea 
individuals were mainly distributed in Hap2, 5, 8, 10, and 
14, which takes 47.92% of the total diarrhea individuals 
(46/96), therefore, those 5 haplotypes may be susceptible to 



Table 6. Effects of different genotypes on piglet diarrhea score 





Exonl 






Exon2 






Exon4 




Genotype 


Number 


LSMiSE 


Genotype 


Number 


LSM+SE 


Genotj'pe 


Number 


LSM+SE 


AA 


34 


0.69+0.18"' 


AA 


30 


1.12+0.19" 


AA 


117 


1.06+0.10' 


BB 


7 


0.83±0.41'"*= 


BB 


59 


1.07+0.14'' 


BB 


9 


0.29+0.35" 


AB 


30 


0.85+0.21*^ 


AB 


89 


0.83+0.11* 


AB 


48 


0.59+0.16" 


BC 


14 


1.77+0.28" 


BC 


25 


0.79+0.21* 


BC 


11 


0.76+0.32"" 


AC 


99 


0.97+0.11" 


AC 


9 


0.22+0.35" 


BD 


2 


0.97+0.74* 


CC 


32 


0.52±0.19' 


CC 


4 


-0.40+0.53" 


CE 
Cf 


23 
6 


1 .05±0.22* 
0.63+0.43"'' 



LSM+SE, Diarrhea score was estimated for each genotype class and shown as means least square means+standard error. 
"■''Values with different superscripts within the same hne differ significantly at p<0.05. 
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Table 7. SLA-DRA gene haplotype and its association withi piglet diarrliea 



Name 


Haplotypes 


Frequency (%) — 


Number 


Diarrhea piglets 

Frequency (%) 


Hapl 


ATACTAAAGG 


8.63 


4 


4.17 


Hap2 


ATACTGAAGG 


4.27 


7 


7.29 


Hap3 


ATCCTGAAGG 


2.60 


1 


1.04 


Hap4 


GTACTAAAGG 


8.89 


3 


3.13 


Hap5 


GTACTGAAGG 


4.13 


8 


8.33 


Hap6 


GTACCAAGGA 


2.08 


4 


4.17 


Hap7 


GTCCTAAAGG 


9.67 


4 


4.17 


Hap8 


GTCCTGAAGG 


8.13 


13 


13.54 


Hap9 


GTCCCAAGGA 


2.64 


0 


0 


Hap 10 


GTCTTAAAGG 


6.25 


10 


10.42 


Hap 11 


GTCTCAAGGA 


3.28 


3 


3.12 


Hap 12 


GTCTCGGAGG 


2.64 


3 


3.12 


Hap 13 


GCACTGAAGG 


6.64 


5 


5.21 


Hapl4 


GCCCTAAAGG 


6.54 


8 


8.33 


Others 




23.61 


23 


23.96 



SLA, swine leukocyte antigen. 



piglet diarrhea; no diarrhea piglets were in Hap9, it may be 
the resistant haplotype. 

DISCUSSION 

Genetic polymorphism of tlie SLA-DRA gene 

Early studies have described limited variation in the 
SLA-DRA locus for the NIH miniature pig (Sachs et al., 
1988; Chardon et al., 1999), but subsequent researches 
revealed that SLA-DRA locus exhibits limited 
polymorphism with 13 alleles representing 3 allele groups, 
compared with other SLA class II genes (Ho et al., 2009; 
Lunney et al., 2009; Ho et al, 2010b). In this study, we 
detected polymorphisms throughout the entire four exon 
regions of SLA-DRA gene firstly, 3, 2, and 6 nucleotide 
mutations were detected in exon 1, exon 2, and exon 4 of 
DRA gene, respectively, which confirmed the result of Ho et 
al. (2009) and Lunney et al. (2009) who showed that the 
SLA-DRA gene is polymorphic. 

Moreover, data on SLA class II genes now predicts that 
the determinants of antigen binding reside in exon 2 
(Chardon et al., 2000). In the SLA-DRA locus, to date only 
two polymorphic sites have been found in the al domain 
(exon 2). One is an A, C, or G (the G has only been found 
in one allele) mutation, which results in a substitution of a 
methionine, leucine, or valine, respectively. The second is a 
C or T which both code for an arginine. The two SNPs 
(c.3093A>C and c.3104C>T) detected in our study are 
accord with the A and C polymorphisms and the C and T 
polymorphism. The only two nucleotide mutations in this 
region may be attributed to the conservative property of 
SLA-DRA gene, these conservative substitutions result from 
the survival of only mutations that retain the the crucial 



function of SLA-DRA gene in regulation of immunity, its 
conservative property is a result of evolutionary selection. 

Only two homozygous genotypes (AA and BB) were 
present in exon 4 of the SLA-DRA gene, alleles C, D, E, and 
F exist in heterozygous forms. There are two factors 
contributing to this: i) there could be a selection effect as 
these alleles have negative effects on individual 
performance, heterozygous individuals have a greater 
advantage over homozygous against a variety of pathogens 
(Sommer et al., 2005), consequently the missing 
homozygous genotypes may have been eliminated in the 
breeding process (Beuzen et al., 2000); ii) there could be a 
sampling effect because the sample sizes were relatively 
small and few homozygous animals with these alleles were 
expected. Thus, this polymorphism mechanism may be a 
necessity for resisting envirormient stress and maintaining 
normal function of SLA-DRA gene. 

Effects of genotypes on piglet diarrhea 

Previously studies have established that piglet mortality 
rate during pre-weaning caused by diarrhea was linked to 
SLA haplotype (Renard and Vaiman, 1989). Our current 
study reported that the diarrhea score of the DQA exon 2 
genotype CC (0.36+0.19) was significantly lower than AB 
(0.98+0.09), BB (0.85+0.77), BC (1.25+0.23) and BE 
(1.19+0.19) (Yang et al., 2013). In the present study, our 
analysis also revealed significant differences in the score of 
piglet diarrhea among different genotypes of SLA-DRA 
gene, this supports the previous observations by Renard and 
Vaiman (1989) and Yang et al. (2013). 

The DRA gene is located at the most centromeric of 
SLA II gene cluster, encodes the a chain with DQA gene. 
The presence of a small number of haplotypes shared by 
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multiple individuals is indicative of linkage disequilibrium 
(LD). Our research group had observed a significant 
association between DQA gene and piglet diarrhea relating 
to nucleotide mutations c.4086A>G and c.4087C>T, 
c.4089A>G and c.4090C>T as linked codons variation in 
allele C of Large White, Landrace and Duroc pigs (Yang et 
al., 2013). According to the results of Yang et al. (2013), 
linkage disequilibrium analysis between the exon 2 loci of 
SLA-DRA-DQA gene revealed five haplotypes of DRA*A- 
DQA*B, DRA*B-DQA*B, DRA*B-DQA*C, DRA*B- 
DQA*D, DRA*C-DQA*B. Absolute value of D' and r^ for 
all SNP loci between DQA and DRA gene showed that of 
the 210 pairwise comparisons, 95 showed significant LD at 
p<0.01 level, indicating the comprising haplotypes of SLA- 
DQA and DRA gene are in complete Unkage disequilibrium 
(ID'I = 1.0000, r^ = 0.7952). In this study, we obtained SNPs 
in SLA-DRA gene and significant genetic correlation with 
piglet diarrhea were found in these 3 observed pig breeds, 
which together may explain the linkage of those two SLA II 
a chain markers reflecting the important function in disease 
resistance. 

Amino acid mutations in coding regions have 
significant impact on the stability of proteins, which largely 
depend on the chemical and physical properties of mutation 
amino acid (molecular weight, volume and hydrophobic), 
mutation position (whether in transmembrane domain), and 
charge properties (Wang and Moult, 2001). The mutations 
in this study were associated with three synonymous 
mutations and seven non-synonymous mutations, in which 
c.l78G>A and c.4167A>G both caused changes of 
uncharged amino acids to positively charged amino acid 
(p.lOGly>Arg, p.206Gln>Arg), and c.211T>C and 
c.4184A>G which resulted in a hydrophilic amino acid 
change to hydrophobic amino acid (p.21Ser>Pro, 
p.212Thr>Ala). Especially, the function of transmembrane 
protein is influenced by single peptide, c.l78G>A that 
occurred in highly hydrophobic areas of the SLA-DRA 
single peptide region, and results in a change of 
hydrophobic residual base to hydrophiUc and therefore may 
lead to the functional loss of this peptide. These are notable 
variations as they not only occurred in an important 
functional region but also caused changes in amino acid 
properties. The constructed haplotypes of all SNPs in the 
three exons and their association with piglet diarrhea 
revealed that haplotypes Hap2 (ATACTGAAGG), Hap5 
(GTACTGAAGG), Hap8 (GTCCTGAAGG), HaplO 
(GTCTTAAAGG) and Hapl4 (GCCCTAAAGG) may be 
the susceptible haplotypes to piglet diarrhea and Hap9 
(GTCCCAAGGA) may be the resistant haplotype. The 
nucleotide substitutions may be linked to variations of each 
other. Thus, we ought to focus on optimizing individuals of 
Hap9 and firstly eUminating individuals of Hap2, 5, 8, 10, 
and 14 in molecular breeding. 



In conclusion, the present study firstly reports five novel 
SNPs within the exon regions of the SLA-DRA gene, which 
will provide more genetic information on the SLA-DRA 
gene variation. Our data revealed that these detected SNPs 
significantly influenced piglet diarrhea, and could be used 
as potential genetic markers for appUcation in pig disease 
resistance breeding. 
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